Biphoton Double-Bell States and Ququarts with "Invisible" Variables 
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We analyze features of mixed biphoton polarization states which arise from pure states of 
polarization-frequency biphoton ququarts after averaging over frequencies of photons. For mixed 
states we find their concurrence C, Schmidt parameter K and degree of polarization P, as well as 
the von Neumann mutual information /. In some simple cases we find also the relative entropy Srei 
and the degree of classical correlations Cd- We show that in mixed states the Schmidt parameter 
does not characterize anymore the degree of entanglement, as it does in pure states. Neverthe- 
less, the Schmidt parameter remains useful even in the case of mixed states because it remains 
directly related to the degree of polarization. We compare results occurring in the cases of full pure 
polarization-frequency states of ququarts, mixed states (averaged over frequencies) and states with 
separated high- and low-frequency parts. Differences between these results can be seen in experi- 
ments with and without a dichroic beam-splitter, as well as with and without frequency filters in 
front of detector. 



I. INTRODUCTION 

As known, a pure state of two particles is entangled 
if its bipartite wave function cannot be factorized, i.e., 
presented in the form of a product of two single-particle 
functions, 



*(xi,X2) X x(a;2), 



(1) 



where xi and X2 arc variables or sets of variables of two 
particles. A natural extension of this definition is given 
by the Schmidt decomposition, or Schmidt theorem [llQ, 
according to which any entangled, unfactorable, bipartite 
wave function can be presented as a sum of factorized 
terms 



^'(a;i,X2) 



A„ (p„{xi) X Xn{x2) 



(2) 



The orthogonal and normalized functions (p„(a;i) and 
Xn{x2) are the Schmidt modes defined as eigenfunc- 
tions of the reduced density matrices pi^\xi,Xi) and 

(2) 

Pr {x2,x'2), which are defined in their turn as par- 
tial traces of the full density matrix /9(a;i, X2; x'j^, ajj) — 

^'(xi, 2:2)4'* (x']^, over X2 (forp^^^) or over 2:2 (forp^^^); 
A„ denotes coinciding eigenvalues of both reduced den- 
sity matrices pf^"* and p?\ n = 1, 2, ... The amount of 
efficiently populated Schmidt modes in the Schmidt de- 
composition ([2]) characterizes the degree of entanglement 
and can be evaluated by the Schmidt parameter 
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For shortening formulations, let us refer the character- 
ization of entanglement by the Schmidt decomposition 
([2]) and the Schmidt parameter K ^ as the Schmidt 
entanglement. Note also that the Schmidt decomposi- 
tion and the definition of the Schmidt parameter (jS]) are 



valid equally for bipartite systems with either discrete or 
continuous variables xi,2. 

It is rather important to analyze the relation of these 
definitions with those arising in the occupation-number 
representation of bipartite states. In the occupation- 
number representation, states are characterized by modes 
1 = 1,2, ...N, and numbers of particles in these modes, 
rii. In the case of bipartite states '^i = 2. Each 
specific distribution of particle occupation numbers in 
modes can be referred to as a configuration and denoted 
as {rii}. The bipartite density matrix in the represen- 
tation of occupation numbers is defined by its matrix 
elements ({fii}|/5|{n^}). On the other hand, each particle 
has its "coordinate", or sets of coordinates, xi and X2, 
correspondingly, for particles 1 and 2. Coordinate ma- 
trix elements of the density matrix p between the coordi- 
nate eigenstates detrmine the coordinate density matrix 
p{xi,X2',x[,X2) = {xi,X2\p\x'i,x'2)- The sets of states 
|a::i,a;2) and are equally complete, and there is the 

following relation between the density matrix in the co- 
ordinate and occupation-number representations 



p{xi,x2;x[,x'2) 



Xi,X2\{ni}) 



x({nJ|/5|K}> (K}|:^'i,4; 



(4) 



The transformation coefficients have sense of the bipar- 
tite coordinate wave functions for given configurations 



{xi,X2\{ni}) 
{{n^\x[,X2) 



*{«i}(a;i,a;2), 
*{„a(^'i'^2)- 



(5) 



In the case of distinguishable particles, the given- 
configuration wave functions ([5|) are always factorized, 
and the Schmidt entanglement can arise only if the 
sum over configurations in Eq. (|4]) contains more than 
one term. This case can be referred to as correspond- 
ing to the configuration entanglement. On the other 



hand, in the case of indistinguishable particles, the given- 
configuration wave functions ([5]) must be either symmet- 
ric or antisymmetric with respect to the variable trans- 
positions (correspondingly, in the cases of two-bozon or 
two-fermion states). Some of such symmetrized given- 
configuration wave functions consist of two terms each, 
and in such cases the Schmidt entanglement can arise 
even if the sum in Eq. Q contains only one term. This 
is the case of the symmetry entanglement. In a gen- 
eral case of multiconfigurational states of indistinguish- 
able particles, the Schmidt entanglement is determined 
by an unseparable superposition of the configurational 
and symmetry types of entanglement. Thus, the Schmidt 
decomposition ([2]) and parameter K ((3|) characterize the 
total amount of entanglement in a pure bipartite state 
with both configurational and symmetry parts of entan- 
glement completely taken into account. 

A simple example of these definitions is the polariza- 
tion state of two photons with different polarizations, 
horizontal (H) and vertical (V) ones, and with the state 
vector and the density-matrix operator given by 

\1h, Iv) = aLal-lO), p = \Ih, Iv){1h, lv\- (6) 

Modes of this state are H and V and this is a single- 
configurational state with uh = 1 and ny = 1. The state 
(|n]) does not have the configurational entanglement. For 
this reason it's often said that the state is not 

entangled at all (see. e.g., Refs. [1, Q)- But the discus- 
sion given above indicates clearly that this is not so be- 
cause photons are indistinguishable bozon particles, and 
their symmetrized wave functions ((5)) 4'i/y((Ti, cr2) and 
^Hv{<^i,<^2) for ttie only given configuration ^ consist 
of two terms each 

Ww(o-i,cr2) = , (7) 

and the same for (cr'j^, cTg)- Here and below we use 
notations ai^2 for polarization variables instead of xi_2 in 
Eqs. ^ and (jH), Note also, that, in contrast to 
distinguishable particles, in the case of indistinguishable 
particles (photons) the variables cti and cti are not asso- 
ciated with any specific particle, they have to be inter- 
preted as the polarization variables of one and another 
of two photons, never known specifically which one of 
them. Also, these variables are not associated with spe- 
cific modes. E.g., as seen well from Eq. ([7]), the variable 
CTi can be related either to the mode H or V. In this 
sense variables and modes are entangled, which means 
that there is no one-to-one correspondence between vari- 
ables and modes. In accordance with Eq. (U) and results 
of Ref. 0, symmetry of the wave function ^'//^(cti, (72) 
^ provides entanglement of the state \1hAv) evaluated 
via the Schmidt decomposition and Schmidt parameter 
K, or concurrence C. This result agrees also with the 
conclusion of Ref. @ . 

However, the same state as given by Eq. ([5]) in the 
basis with horizontal and vertical axes Ox,Oy (0°,90°), 
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in the basis turned for 45° takes the form 

|245°,0i35o) — 10450,21350) (ajgo — al35o)|0) 

71 = 2 ■ 

Two modes correspond now to photon polarizations 
along the 45°- and 135°-axcs, and, obviously, the state 
([8]) is a two-configurational one. The configurations are 
(71450 = 2, ni35o = 0) and (71450 = 0, 71.1350 = 2). The sum 
over these two configurations in Eq. (|4]) provides now 
the configurational entanglement. In the same time, the 
wave functions ([5]) of the configurations (2450,61350) and 
(0450,21350) are factorized, and in this case there is no 
symmetry entanglement. Hence, in the basis (45°, 135°) 
the Schmidt entanglement arises only owing to configu- 
rational entanglement, but both the degree of entangle- 
ment evaluated by the Schmidt parameter K, and the 
Schmidt modes remain the same as found in the ba- 
sis (0°,90°). Thus, with changing bases, symmetry and 
configurational entanglement of pure polarization bipho- 
ton states are transformed into each other. In extreme 
cases of the bases (0°,90°) and (45°, 135°) one of these 
two types of entanglement completely disappears. In all 
intermediate cases (a, a + 90°) both types of entangle- 
ment exist and have to be taken into account together. 
The Schmidt decomposition, Schmidt modes, and the 
Schmidt parameter K characterize the total entangle- 
ment with both the symmetry and configurational con- 
tributions taken into account, and this makes the total 
entanglement basis-independent {K{a) ~ const. = 2). 
The total entanglement is the entanglement of a state as 
a whole, independent of conditions of its measurement or 
methods of its theoretical description. 

One general remark more concerns spacial separation 
of indistinguishable particles, which is considered some- 
times as partially diminishing their indistinguishability, 
making particles as if somewhat distinguishable. In the 
frame of such approach indistinguishable particles are as- 
sumed to get some features distinguishing them, e.g., as 
photons propagating to the left and to the right, with 
left- and right-propagating photons treated as usual dis- 
tinguishable particles. Then, in accordance with Eq. 
only configurational entanglement is taken into account. 
In fact, such procedures are insufficient for finding total 
entanglement in systems of indistinguishable particles. 
Actually, any kind of a spacial separation creates a new 
degree of freedom for particles, but particles themselves 
remain indistinguishable. Bipartite states with several 
degrees of freedom (e.g., polarizations and propagation 
angles of photons or polarization and frequencies) are 
more complicated, and have a higher dimensionality than 
states with one degree of freedom (purely polarization 
biphoton states). 

Polarization-frequency or polarization-angle ququarts 
are considered usually as two-qubit states, in which fre- 
quencies (or propagation angles) of photons are assumed 
to be strictly related to their variables, and symmetry 
of the wave function is not taken into account Q. In 
terms of discussed above configurational and symmetry 



3 



entanglement, the two-qubit model of biphoton ququarts 
describes only the configurational entanglement and ig- 
nores the symmetry entanglement inevitably occurring 
in the case of indistinguishable particles. Though such 
approach can be reasonable in some cases, as said above, 
it is insufficient for evaluation of the total entanglement 
of biphoton ququarts. The approach appropriate for 
evaluation of the total entanglement was suggested in 
our work Q. In this approach biphoton polarization- 
frequency ququarts are considered as two-qudit states 
with the dimensionality of the one-photon Hilbert space 
d = 4. Their wave function is symmetric with respect 
to the transposition of photons, and their entanglement 
is a combination of the polarization and frequency en- 
tanglement, as well combination of the configura- 
tion and symmetry entanglement. Features of frequency- 
polarization two-qudit biphoton ququarts are reminded 
briefly below in section II. In section III we discuss a rep- 
resentation of a general ququart's wave function in the 
form of sums of double-Bell states, where "double-Bell" 
means a product of, e.g., polarization and frequency Bell 
states. In sections IV and V wc describe mixed biphoton 
states with "invisible" variables which arise from bipho- 
ton ququarts after averaging over either frequency or po- 
larization variables. In experiment this corresponds to 
the situation when experimenters use detectors unsensi- 
tive either to values of frequencies of photons or to their 
polarizations, which makes these variable effectively in- 
visible. Features of the arising two-frequency or two- 
polarization mixed states are rather unusual. Wc find 
explicitly their Schmidt parameter K, concurrence C, 
von Neumann mutual information / and, in some simple 
cases, their relative entropy Crei and the degree of clas- 
sical correlations Cd- For biphoton mixed polarization 
states we find a universal relation between their degree of 
polarization and the Schmidt parameter K. On the other 
hand, we find that in the case of mixed polarization states 
the Schmidt parameter K does not characterize anymore 
their degree of entanglement as it does in the case pure 
bipartite states. At last, in section VI wc discuss in 
more details a relationship between features of polariza- 
tion biphoton ququarts with invisible frequency variables 
and of two-qubit ququarts with asymmetric wave func- 
tions and photon frequencies not considered as variables 
at all. We discuss experimental schemes in which bipho- 
ton polarization-frequency ququarts can display either 
features of full-dimensionality two-dudit states, or of a 
two-qubit model, or of mixed polarization states aver- 
aged over frequencies. 



nonlinear birefringent crystal and/or additional manip- 
ulations with photon polarizations after crystals. Given 
frequencies of photons in SPDC pairs (oj/j and uji) can 
be obtained either with the help of frequency filters or 
in the approximation of a highly monochromatic pump 
and a very long crystal. In such schemes, polarizations 
and frequencies of photons are two independent degrees 
of freedom, and the corresponding variables ai and uji 
{i — 1,2) can take two values each for each photon, 
ai = H or V and uJi = LUh or uji. The only restric- 
tion is that frequencies of two photons are always dif- 
ferent, LUi ^ UJ2- One-photon polarization- frequency ba- 
sis modes are iJZ, Vh, andT^L The corresponding 
state- vectors and wave functions are given by 



state 
vector 



wave function 



pol 



pol 



pol 



freq 



freq 



freq 



freq 




(9) 



Here the superscripts pol and freq refer the polarization 
and frequency degrees of freedom. The upper and lower 
lines in the two-line columns correspond to the horizonal 
polarization and high frequency, whereas the lower lines 
- to the vertical polarization and low frequency. Alto- 
gether, four one-photon states Q form a qudit with the 
dimensionality of the one-photon Hilbert space d = 4. 

Under the assumption that frequencies of two photons 
in the SPDC pairs are always different one can construct 
from four one-photon states © only four independent 
basis biphoton states with the state vectors 



I* 



HV/ ~ "-Hh^Vll 



(4) \ _ „t J 
HV/ ~ " 



-|0), 



(10) 



and corresponding to them wave functions [3] 



II. BIPHOTON QUQUARTS 

Biphoton ququarts can be produced, e.g., in collinear 
non-degenerate SPDC processes [two possible frequen- 
cies, uJh (high, h) or (low, I) and one of two polariza- 
tions for each photon, horizonal (H) and vertical {V)]. To 
get ququarts of a general form, on can use more than one 




(11) 



4 



as 



(4) 
HV 



1 

71 



:)r«0r«(;)r«0 



)frec 




(12) 





free 



(14) 



where the indices 1 and 2 numerate variables of indistin- 
guishable photons (7i , wi and (T2, a;2 (not seen explicitly in 
the matrix form of writing the wave functions). All basis 
wave functions (jll[) - p4p are symmetric with respect to 
the transposition of particle variables 1 ^ 2, as it should 
be for any two-bozon states. In a general form, the wave 
function of a ququart is defined as a superposition of four 
basis wave function ([TT |) - (|14p 



,(4) 
HV 



,(4) 
VH 



,(4) 
VV 



(15) 



freq 



^ HH — 



(4) 
HV 



^ HV — — 




^VV ^ 



/0\ 








(?) ^ 









/req 



(18) 
(19) 



Biphoton basis states written in the form of Eqs. (fTTj)- 
(|14l) are convenient for reduction of the density matrix 
over one of two pairs of photon variables, (ui, wi) or 
(ct2, 1^2), whereas the form (|16p - p^ is more convenient 
for reduction in one of two degrees of freedom, polariza- 
tion or frequency, i.e. for transition to mixed states. 

As evident from Eqs. (ITT|)-([14| and (HHl-dTll), bipho- 
ton ququarts describe two-qudit states of dimensionality 
D = cP = 16, and they can be considered as two-qudit 
states (with d = 4) Q. This understanding differs from 
a widely spread opinion that biphoton ququarts are two- 
qubit states. The two-qubit model is correct only for 
ququarts of two distinguishable particles with only one 
degree of freedom taken into account. In contrast to this, 
in the case of biphoton ququarts we have two equally 
important degrees of freedom (e.g., polarization and fre- 
quency) and, owing to obligatory symmetry of two-bozon 
states, entanglement in both degrees of freedom is possi- 
ble and has to be taken into account. This makes features 
of pure-state biphoton ququarts rather significantly dif- 
ferent from those of the two-qubit ququarts . 

Note also that instead of frequencies LOh and lui one can 
consider photons with coinciding frequencies but differ- 
ent directions of propagation (in the non-coUinear SPDC 
scheme) characterized by angles Oieft and Bright- In this 
case two degrees of freedom are polarizations and angles, 
determining directions of propagation of photons, and 
two sets of variables are (cti,6'i) and (cr2:^2)- Except 
different notations, the physics of polarization-frequency 
and polarization-angle biphoton ququarts is identical. 



with arbitrary coefficients Ci_2.3,4 obeying the normaliza- 



tion condition |Ci| 



IC2 



|C4K = 1. 



Note that in Eqs. ([TT|) - ([T4l) the two-row columns are 
combined into four-row ones corresponding to given par- 
ticle variables but mixed up polarization and frequency 
degrees of freedom. Alternatively, it's possible to com- 
bine parts of wave functions with given degrees of free- 
dom but mixed up variables to rewrite Eqs. (|lip - P^ 



III. DOUBLE-BELL STATES 

Instead of four basis wave functions of ququarts pT|) - 
T4| one can use their superpositions: 



(4) 



1 
V2 



$(4) ^^(4) 
^ HH ^VV 



(4) _ 1 



(4) 
HV 



± ^ 



(4) 
VH 



(20) 



5 



The sum of these basis functions with arbitrary coeffi- 
cients gives a general expression for the wave function of 
a biphoton quaquart equivalent to that of Eq. ()15p 



(4) 



(4) 



(4) 



(4) 



(21) 



where C± = (Ci ± C4) /V2 and B± = (C2 ± C3) /V2. 

Below we will use also a mixture of basis wave functions 
(Ulll-dll and (1201): 

{^'Sh. *^'^}- (22) 

Expansion of a general ququart's wave function in the 
basis functions of Eq. is given by 



= + + C^^^y\ + (23) 

with 



|Ci| 



IS4 



(24) 



The wave functions and $± ^ of Eq. ([201) are not 
Bell states because they have a higher dimensionality 
{p = 16) than the 4-dimcnsional two-qubit true Bell 
states 



{Bell) 



(Bell) 



1 r 

V2 



1 

V2 



± 



± 




But, as can be easily checked, the functions 5'^. and $j_ 
can be presented as direct products of polarization and 
frequency Bell states, and these products can be referred 
to as double-Bell states: 



(4) 



(4) 
± 



(Bell)pol 
± ' 

(Bell) pol 
± 



(Bell) freq 



(Bell) freq 



(27) 



With these equalities taken into account, the general 
wave function of a biphoton ququart (j23|) can be rewrit- 
ten as 



Bell) freq 



(Bell) pol 



(Bell) freq 



where 



(28) 

^i3)poi ^ Ci^P^H + B+^^'^^P"' + C^^Py^l (29) 



is a general wave function of a polarization biphoton 
qutrit 01 ; 



pol 
HH 



pol 
1 



pol 



pol _ (0 
VV ' 



pol 



pol 



(30) 



are purely polarization states of two photons with coin- 
cidingpolarizations. For comparison with notations of 
Ref. Q, and C4 in Eq. (^9)) correspond to C2 and 
C3 for qutrits in 

The density matrix of a ququart (P5| is given by 

„(4) _ 



^(3)po; ^ ^iBell)freq _^ g_^{Bell)pol ^ ^(B ell) freq 



(Be;/) freq 



B^^ 



(Bell) pol 



{Bell) freq 



(31) 



As mentioned above, the dimensionality of the density 
matrix p'"*^ is 16 x 16. But the number of ququart's basis 
states in all equations (fT0)) - (|28)) is only 4 rather than 16. 
This means in fact, that 12 missing basis states come in 
the superpositions forming the biphoton ququarts with 
zero coefficients. The excluded states are states with a 
wrong symmetry and states with coinciding frequencies 
of photons. Owing to these omissions the 16 x 16 den- 
sity matrix has many zero elements. Positions of these 
zero elements depend on a choice of a basis used for 
calculation of the density matrix. The density matrix 
has the simplest form if it is calculated in the basis of 
ququart's basis states, e.g., in the basis of states (|22p . 

1*^1^, ^yv-} plus 12 other unspecified ba- 

sis states coming with zero coefficients in the ququart's 
wave function (|23p . With such choice of the basis for 
calculation of the density matrix, its 4 nonzero lines and 
columns can be concentrated in the upper left corner. 
Then, the nonzero 4x4 part of the density matrix ([?T|) 
can be written as 



IV. 



|Ci|2 Cis; C1C4* CiB*_ 

B+Cl 15+ 12 B+Cl B+Bl 

040]^ 04-D_]_ |04 1 04-D_ 

\B^Cl B^BX B^C'l \B^\^J 



BIPHOTON QUQUARTS WITH 
"INVISIBLE" VARIABLES 



(32) 



For measuring parameters of biphoton ququarts one 
has to use detectors (photon counters) provided with 
polarizers and frequency filters to distinguish photons 
in all four polarization- frequency modes (jOl) 01 • How- 
ever, it's interesting also what happens if one uses only 
part of these devices, either only polarizers or only fre- 
quency filters. In these two cases either frequency or 
polarization variables become "invisible". Mathemat- 
ically, such states are characterized by density matri- 
ces pP°'' = Trfreqp Or pf^'^'i = Trpoip, where p is the 
general density matrix of Eqs. (pij) . (|5^ . and Trf^eq, 
Tvpoi denote partial traces over frequency or polarization 
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variables. Features of such states with "invisible" vari- 
ables are significantly different from those of both original 
biphoton ququarts and purely polarization or frequency 
biphoton states. If original biphoton ququarts are pure 
states, biphoton ququarts with invisible variable are, typ- 
ically, mixed. Besides, if density matrix of the original 
biphoton ququart is 16-dimensional, the matrices 
and pf'^'^'i have the dimensionality D = A and describe 
in a general case two-qubit mixed states. Features of 
such states are rather interesting and they are described 
below. As the first example, let us consider the two- 
frequency states with "invisible" polarization variables. 



Two-frequency state with "invisible" polarization 
variables 

The density matrix pf^^'J = Tvpoip is easily found from 
Eq. ([3T|) to be given by 



jfreq ^ (1- 



Bell) freq 



(Bell) freq 
+ 




nat. basis 



(33) 
(34) 

(35) 

ell basis 

Eq. (|34p corresponds to writing the frequency den- 
sity matrix in the natural basis < (0)70)1(1 



whereas Eq. ([35]) - to its presentation in the basis of Bell 

states ^(Bell)freq ^{B ell) freq ^(Bell) freq^ 

with zero lines and columns dropped. 

As it's evident from Eqs. and (|55|) . eigenvalues 

and entropy of the matrix pf^'^'i are equal to 



Ai = A2 = 0, A3 = 1 - A4 = |B_ 



(36) 



and 



h loga A 



--2|i?_|^log2|i3-| 



(l-|B_P)log2(l-|i?-n. 



(37) 



0.5 



S 



0.2 0.4 0.6 2"" 1 

Figure 1: Entropy I I37I I of the two-frequency mixed state with the 
density matrix pi^'^i | |34ll . II38I I as a function of \B—\ 



To calculate the reduced density matrix of the mixed 
two-frequency state, we have to present the full matrix 
pf^'^t ((55)) in the form of a sum of products of 2 x 2 
matrices 



P 

1 




freq 



E 

± 





1±(1-2|S_|2) 



± 



1/,^ vo 0/1 



4 

1 





1 



± 





1 



1 





1 



1 




(38) 



By taking the trace of this matrix, e.g., with respect to 
the frequency variable of a photon 2, we easily find that 
the reduced density matrix is given by 



„J'req 



1 fl 



2 VO 1 



(39) 



Eigenvalues of this matrix are are A^'^'' = X^"^ = i, 
which corresponds to the Schmidt parameter K^^'^'' = 
[A^'^''^ -I- Aj*^-*^]"^ = 2 and entropy of the reduced state 
S{p(^'^'^) = -2 X i X log2 5 = 1- By combining the last 
result with that of Eq. p7l) . we find the von- Neumann 
mutual information of the mixed frequency state charac- 
terized by the density matrix of Eqs. (p3|) . (|34|) 



lipf'-"'^) = 2S{pl''"') - Sip^''"'^) 



Sipf'''"') 



= 2+ ii?_piog2 + (1 - ii?_niog2(i - 



(40) 



The entropy S{p^^'^'^) is shown in Fig. [T] as a func- 
tion of It's seen to turn zero only at B- ^ 
and = 1. In these two cases the density matrix 
pfreq ^escribes pure states: at i3_ = Eq. ([55]) yields 

freq ^ ^(Bell)freq ^ 



(Bell) freq 
+ 



pfreq _ ^(Bell) freq ^ ^^{Bell) freq^ 



and at |£?_| = 1 



In all other cases, 



\B- 1^0,1, the state characterized by the density matrix 



P 



freq 



dSSD-dSni) is mixed. 



in agreement with the results of Ref. Q. The concur- 
rence of the mixed state under consideration can be easily 
found directly from the first expression of Eq. with 
the help of the Wootters' procedure @, and it's given by 



fjfreq 



Il-2|S_| 



(41) 



In addition to the concurrence, there are several other 
measures of entanglement valid for mixed states. One 
of them is the relative entropy suggested by Vedral et al 
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[8|, [lOl, which is defined as the "distance" between 
the density matrix p and the closest unentangled density 
matrix cr: Srei{p) = Tr[p{\og2 p — logj u)]. For the state 
of the form ([55)1 the relative entropy is known Q and in 
a slightly modified form it can be presented as ((33| 



1 + C.f-""! 



I0g2(l + C^"') 



1 _ (jfreq 



l0g2(l-C^'-'=«). 



(42) 



At last, one possibility of quantifying the amount of 
classical correlations is related to the definition of works 
0, [13, where the degree of classical correlations was 
defined as the difference between the von Neumann mu- 
tual information and the relative entropy, for the state 
([33)1 given by 



(43) 



again, in agreement with the conclusion of Ref. [8|. All 
characteristics of correlations in the mixed two-frequency 
state are shown in Fig. [5] as functions of the parameter 
\B-\. This picture shows that: 1) The relative entropy is 




0.5 



0.2 



0.4 



0.6 2" 



\B\ 1 



Figure 2: Parameters, characterizing degrees of correlations in the 
two-frequency mixed state with the density matrix pf'^'^i ||34l l as 
functions of | -B— | . 

always less than the concurrence except the cases of pure 
states = or 1) when they are equal; hence, in the 
case of mixed states the relative entropy better character- 
izes the level of quantum correlations than concurrence. 
2) In the case of pure states entropic characteristics of 
the levels of classical and quantum correlations are equal, 
and each of them equals a half of the von Neumann mu- 
tual information, = C^'~'=« ^ Sr ^ ^Hp^'^'"') in the 
cases B- = and \B-\ = 1. 3) The Schmidt parame- 
ter K is a non-cntropic characteristic which can be used 
for evaluation of the levels of both classical and quan- 
tum correlations only in the case of pure states; in the 
case of mixed states K{\B-\) behaves as Cci{\B-\) (they 
are both constant) and, hence, for mixed states, K can 
be considered only as a non-entropic measure of classi- 
cal correlations; relation between the Schmidt parameter 
and concurrence C = -^2(1 — K~^) occurs only in the 
case of pure states and is inapplicable for mixed states. 



V. MIXED POLARIZATION BIPHOTON 
STATES WITH "INVISIBLE" FREQUENCY 
VARIABLES 



If we decide to choose the frequency variables of pho- 
tons as "invisible" ones, we have to take a trace of the 
general density matrix p^'*^ (|3ip with respect to "invis- 
ible" frequency variables to get the density matrix of a 
mixed two-photon polarization state 



(44) 



Reduced density matrix, eigenvalues, Schmidt 
paramerter, and entropy 



For finding further the reduced density matrix of the 
mixed polarization biophoton state, let us write down 
first the density matrix p'P°^ in the form of products of 
single-particle matrices 



\Ci\ 

1 




1 




1 


1 



1 




|i?4 



1 

71 



CiB* 



/O 1 

+ U 



1 

'71 



B+Cl ( Q J 



2 


1 



+ CIB+ 
B^Ci 



2 


1 



1 


1 



1 



1 




2J 





1 





'I 




® 1 



















(45) 



In this expression only a half of terms is written explic- 
itly whereas the other part is indicated by ... and, actu- 
ally, dropped as unimportant for calculation of the den- 
sity matrix reduced over the variable 2, because in these 
terms all matrices (* are off-diagonal. Thus, the re- 
duced density matrix determined by Eq. (|45p is given 
by 



\Ci\ 



|g+r+|g-r 



72 



C 1 C4 



where 



\Ci\ 



\B\ 



z 

1-x 



(46) 



CiBl + B+Ct 
V2 



(47) 
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Eigenvalues of the matrix p6|) are equal to 

,r)pol _ 1± V4|zP + (2x-l)2 
A-l- — T 



l±J{l-\B_\^f-\2C,C,-Bl\^ 



The equality 

4|zp + {2x - I f = (1 - - I2C1C4 - B 



(48) 



2 |2 



is proved with a simple algebra [with the use of the nor- 
malization condition (|24p] and it shows, in particular, 
that at any values of the ququart's parameters 



1- > \2CiCi~Bl\ 



(49) 



The Schmidt parameter, determined by the eigenvalues 
is given by 



1 



1 + (1- |S_|2)2_ |2CiC4-S2|2 



(50) 



Entropy of the reduced mixed polarization state is given 
by 



5(pr) = l-E 



l±./(l-|i3_|2)2_|2CiC4-S2|2 



log2 1 ± V (1 - IS- P)- I2C1C4 -Bl\^\. (51) 



B. Degree of polarization and the Schmidt 
parameter 

The reduced density matrix (|46p can be used also for 
finding the one-photon polarization vector, or vector of 
Stokes parameters per photon 

^=Tr {pP°^a) = {2Re{z), -2Im{z),2x - 1}, (52) 

where a is the vector of Pauli matrices. In terms of Stokes 
parameters the degree of polarization is determined as 



P 



where S is the biphoton polarization 

vector [17|. By comparing Eqs. ([50)) and ([5^ we find 
the following general relation between the degree of po- 
larization and the Schmidt parameter 



P^ + 2 1 - {KP°^y 



= 1. 



(53) 



This is the analog of the relation + ~ 1 found 
in Ref for pure states of biphoton qutrits, for which 
2{K~^ — 1) = (72. In the case of mixed states we consider 
here, their concurrence C^°' (found below) is not related 



anymore with the Schmidt parameter Kp°K But the re- 
lation (|53p between the degree of polarization and the 
Schmidt parameter holds good. This shows once again 
that for mixed states the Schmidt parameter K char- 
acterizes in some way the degree of occurring classical 
rather than quantum correlations, because the degree of 
polarization is a purely classical quantity. 

Note also that the degree of polarization in mixed po- 
larization states coincides exactly with the degree of po- 
larization of the original two-qudit quqart as a whole, 
Pmixed = P^^'' = \Tr (^p^^^a)\. This is clear because to 
calculate P*-"*^ we have to take the same traces of the 
ququart's density matrix p^^^ (PT|) . ([5^ which reduce it 
first to dUl), then to pP°' (gB]), and in this way to 
the definition ([5^ of the vector of Stokes parameters. 
This procedure reduces the definition of P*-^' to that of 

Pmixed- 



Full density matrix: eigenvalues, entropy, and 
von Neumann mutual information 



Returning to the general density matrix of 
a mixed polarization state (|44p . explicitly it 
takes the simplest form in the basis of states 



dSni) and (1201): 



^pol 



pi3)pol 



IP- I 



/|Ci|2 CiP; CiCi 

B+Cl |P+|2 B+Cl 

C4C1 Cj^B*^ IC4P 

V 



\ 




|P_|2/ 



(54) 



where is the qutrit's coherence matrix (l5l - [l7| . 0. 

By comparing with the full ququart's density matrix of 
Eq. ([5^ we see that averaging over frequency variables 
obliterates terms linear in P_ and PI. The expression 
on the right-hand side of Eq. ([54)) was obtained also in 
the paper [Tsl , though not in the context of mixed polar- 
ization states averaged over the frequency distribution. 

Eigenvalues of the density matrix pP"^ ([5l)) can be eas- 
ily found and, interesting enough, they appear to coin- 
cide with eigenvalues of the two-frequency mixed state 
considered above ([55)) . Ai^2 = 0, A3 = |P-P, and 
A4 = 1 - |P_|2. Hence, the entropy S{pP°^) IS given 
by the same equation as S{p^^'^'^) ([37)) and in its depen- 
dence on |P-| is characterized by the same curve as in 
Fig. [TJ Again, the polarization state with invisible fre- 
quency variables is pure only in two cases, P_ = and 
|P_| = 1, and otherwise it is mixed. 

A special note should be done about the case |P_| = 1 
(and Ci = P+ = C4 = 0); when the polarization 
density matrix of Eq. ([44)) is reduced to pP°'|b_=i = 



{BeU)pol 



(* 



{BeU)pol 



This density matrix de- 
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scribes a pure polarization state with the antisymmet- 
ric wave function ^^°'|b_=i = \^iBeii)poi ^ r^j^^ result 
can seem to be in contradiction with the general re- 
quirement of the Boze-Einstein statistics according to 
which all biphoton wave functions are obliged to be 
symmetric. But, actually, there is no contradiction be- 
cause in the case under consideration the polarization 
wave function is only one part of a more 

general polarization-frequency biphoton wave function 

^BPh\B_=l = ^ff"'"^^"' ® ^(BeU)fre,^ ^j^.^j^ SymmCt- 

ric. In a more general formulation, to produce antisym- 
metric pure polarization Bell state vp^^'^^'P"' one has to 
produce, in fact, a more general symmetric state with 
at least one degree of freedom additional to polarization 
such that (i) the total wave function is factorized, i.e, 
is given by a product of parts depending on polarization 
and other degree(s) of freedom and (m) both parts are an- 
tisymmetric. Then by ignoring the additional degree(s) 
of freedom one can observe the antisymmetric pure po- 
larization state \]>(^'^'')J'°'_ In the case of the polarization- 
frequency biphoton ququarts, the role of such additional 
degree of freedom is played by photon frequencies. 

The von Neumann mutual information of the mixed 
polarization state is determined by the usual formula 



(55) 



with SipP"') and 5(pP°') given by Eqs. 1^ and §7\ 



where 



p(3)^(3) ^ |^(3)<)^^(3)|^(3)<)^^(3) 
= -(2C7i*C4* -B;2)|*(3)^^$(3)| 



--B+C4 B, 



B+Ci 



(58) 



— C4 CiiB^ C4C1 



Eigenvalues of the matrix pP°^pP°^ are given by j"* plus 
three eigenvalues of the 3D matrix p^^^p'-^\ Eigenvalues 
of the last matrix can be found explicitly to yield 



(59) 

Square roots of these values are just the Wootters eigen- 
numbers: 



^(2 



A3) 



(4) 



^(2) 



O7 



IB- 1 



(60) 

Concurrence of the mixed polarization state is given by 



the difference between 



and the sum of all other 



A 



w 



CP°^ = \ \2CiCi~Bl\ - |B_ 



(61) 



D. Concurrence of the mixed polarization states 

Concurrence of the mixed polarization state is deter- 
mined by means of the standard Wootters procedure Q . 
It's best starting from for the density matrix pP°'' written 
in the form of Eqs. (|54)) . Note also that the qutrit's co- 
herence matrix in the first of these equations can be pre- 
sented as = |^r(3)^)^^(3)| |^(3)^) ^ ^^^^ rpj^g 

next step consists in constructing the spin-flipped and 
complex conjugated matrix pP°K Spin flipping means 
changing polarizations H ^ V with multiplication of the 
wave functions by i for H- and —i for l/-polarization. 
This gives the following transformation rules for the 
ququart's coefficients Ci — > — C4, i?+ 5+, C4 — > — Ci, 
B- — >■ —B-. As a result we get the following expression 
for the matrices p and p^^^ 

P |B_|2j' 

p(3) = |$(3))($(3)|, |$(3)) = (^B+y (56) 

With the help of these expressions, the product of matri- 
ces pP°^ and pP°^ takes the form 



E. Example 1 

Ci = C4 = 0. This is the case when the original 
ququart's wave function has the form of a sum of two 
double-Bell states 

^"^^ = 5_^^(Sen)poi ^ ^{Bell) freq 

with l-B+P + = 1. The mixed-state polarization 

density matrix pP°^ (j44p corresponding to this wave func- 
tion is reduced in this case to the form identical to that 
of the mixed-state frequency density matrix pf^^'^ ([55)) : 

Consequently, the states are characterized by the 
same correlation parameter as shown in Fig. [21 In par- 
ticular, in this case RP"'^ = 2 and C^"' = |l-2|B„|2|. 
Besides, owing to the relation between the degree of po- 
larization and the Schmidt parameter ([55[) . P = 0, i.e. in 
this case all states (j63p are unpolarized, independently of 
values of the complex constants i?+ and i?_ . 
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F. Example 2 

Bj^ = 0. In this case the density matrix p^"' takes 
the form 



(Bell) pol 



(64) 



where ^'^j^ and are given by Eqs. ((30|) . Expres- 

sions for the nonzero eigenvalues of the reduced density 
matrix A^'^"', Schmidt parameter K'''°\ degree of polar- 
ization P, and concurrence (7^°' follow from general Eqs. 

dMl), O, dSSl), (ED) 



A 



{r)pol 



i±||CiP-|C4H 



2 ' i + (|Ci|2-|C4Pr 

P=\\C\\^- \Ci\^l CP"' = I (|Ci| + IC4I)' - 1| (65) 



with the normalization condition |Cip -I- IC4P = 1 — 
These expressions are two-parametric, and for 
their simpler analysis let us consider two one-parametric 
cases. 



1. B+^Ci^O, \Cif = l~\B-f 

This case corresponds to the following original 
ququart's wave function 



^{Bell)freq^ (66) 



or in the traditional form p5)) to coefficients C2 
B^/V2 and C3 = -B^/V2: 



(Bell) pol 



(4) 



(67) 



Under these conditions the density matrix of the mixed 
polarization state (pi)) . takes the form 



p^'°' = (i-iB_nvi/^^°i,®(*g'^)^ 

and Eqs. (|65|) are reduced to 



(Bell) pol 



(68) 



, (r)poi 



1 ± 1 - |B_ 



2-2|S_|2 + |S_ 



P=1-\B\' = 1- CP°\ CP°' = |S_ 



(69) 



Note that the derived relation between the degree of po- 
larization and concurrence P = 1 — CP"'' is specific for 
the case under consideration, B^ = C4 = 0. But this 
result is rather important because it shows that the de- 
gree of entanglement can be found in this case via direct 



experimental measurements of the degree of polarization 
of the state ((55)) . 

At B-i- = C4 = the entropy of the reduced state ([68|l 
equals to 



1-E 



S{pP"') = 

1 ± ll 



5]AW^'"log2(Ai^)^'") 



± 

IB- PI 



■log2(l±|l-|B-| 



(70) 



± 



The functions KP"'{\B_\), CP"'{\B^\) and P(|B_|) are 
shown in Fig. |3] together with the von Neumann mutual 
information /(|S_|) = 2S{pP"') - S{pP"') determined by 
Eqs. dSZl) and The left and right borders of this 

Figure correspond to pure two-qubit polarization states, 

1. The state 



at B 



and at m 




0.2 



0.4 



0.6 



IB- 



Figure 3: The Schmidt parameter K^°'- , concurrence C*'"', von 
Neumann mutual information 1^°^ and the degree of polarization 
P as functions of the mixing parameter _B_ for the state 1 I68II . 

^ffff is maximally polarized and has no correlations at 
all, either quantum or classical ones. The state vp^*^"^^"' 
is unpolarized and has maximal possible degrees of classi- 
cal and quantum correlations, C = 1 and Cd = I—C = 1. 
At intermediate values of |i?-|, 1 > |-B-| > the po- 
larization state (pS)) . ([5^ is mixed and possesses both 
quantum and classical correlations. At |i?_| 7^ 0, 1 the 
degree of classical correlations Cci is somewhat larger 
than /(|B_|) — C{Bi) because the curve of the relative 
entropy Srei (not shown here) inevitably is located some- 
what lower than C(|i3_|). 



8. B4 



0, |Ci 



1-|S_|2 



In this case the general density matrix of the mixed 
polarization state (|44|) has the form 



P' 



pol 



l-|i?„ 



^ ^HH + ^ ^VV 

(^BelLpoiy ^ (71) 



\2^Bell,pol 
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where ipi^4 are phases of the coefRcients Ci^. In terms 
of entanglement and correlations, the state ([TT]) is equiv- 
alent to the two-frequency mixed state (|33| and to the 
above considered case Ci = C4 = 0, ^ Eqs. 
(l65l) are reduced to 



. (r)po/ 



1 



^poi 2, P = 0, CP°' = |l - 2|B_|2|, 



(72) 



and all correlation parameters as functions of \B-\ are 
characterized by the same curves as in Fig. (3] 



the density matrix pr arising from the representation (|73p 
the off-diagonal terms disappear and pr takes the form 

(75) 



where 



^ ( |Ci|2-f CiCl+C2Cl\ 



(76) 



VI. COMPARISON WITH THE TWO-QUBIT 
MODEL OF BIPHOTON QUQUARTS 

The consideration given above assumes experimental 
investigation of ququarts with devices nonselective either 
with respect to frequencies or polarizations of photons. 
In the case of spectral mixed states the detectors and 
beam splitters to be used are assumed to be spectrally 
non-selective. On the other hand, there is an alternative 
method of measurements. It's possible to split the bipho- 
ton beam for two channels by a dichroic beam splitter 
to get two beams with only high-frequency or only low- 
frequency photons. After this it's possible to investigate 
spectral correlations only in one of these two channels 
independently of the other one. 

A theoretical picture corresponding to this procedure 
is based on the observation that the general wave func- 
tion of a biphoton ququart (jlSp can be written in the 
"coordinate" representation in the form 



and 



'f'^'*^ (0-1,^1; 0-2, W2) 



where and 



(73) 



2 ^f, oiiu ^ 2 gb ^'^'^ ^^"^ often used two-qubit purely 
polarization ququart 's wave functions with a broken sym- 
metry 



Cl 



pol / \ pol / \ pol 



pol 





,(0 



pol 
2 





pol 



pol 
2 



(74) 



and the same for ^'2^^ with the replacement C2 ^ C3. 

Such wave functions, asymmetric with respect to the 
variable transposition 1^2, cannot exist by themselves, 
without multiplication by frequency parts in the total 
wave function (|73p . But it's possible to think that the use 
of a dichroic beam splitter in experiments can give sense 
even for such "non-existing" wave functions. Indeed, in 



^*(')^(ai,a2)*(';;(a;,a2) 



Cap CiC^+CsCl 



(77) 



A half sum of these partial reduced density matrices give 
exactly the above derived reduced density matrix of the 
mixed polarization state (j46|) 



pol 



Pr, h ^ PrJ 



(78) 



Separation of high- and low-frequency parts in the re- 
duced density matrix pi^'' (|75|) gives rise to a possibility of 
interpreting each of two terms in this equation as describ- 
ing correlations in each of two channels after the dichroic 
beam splitter, independently of another. Note however 
that separation for low and high frequencies arises only 
in the reduced but not in the full density matrix. The 
latter describes all biphoton beam as a whole, either split 
in a dichroic beam splitter or not. 

Eigenvalues of the density matrices pr,h and pr,i (|76p . 
(|77p are easily found to give the well known concurrence 
and Schmidt parameter of a model two-qubit ququart @ 

C2gb = 2IC1C4 - C2C3I - \2CiCi -Bl + B''_\, (79) 

K2qb = ^ L = : ^, (80) 

whereas in mixed polarization states these parameters 
are given by Eqs. (|6ip and ([50]) : 



CP°' = \ \2CiC4-B^ 



\B. 



1 + (1- |B_|2)2_ |2CiC4-S2|2 



(81) 



(82) 



Differences between Eqs. ^ and (EU, dHU) re- 

flect a general difference between the mixed polarization 
state of the ququart as a whole and its part in sepa- 
rate high- or low-frequency channel arising in the reduced 
density matrix. Probably, the simplest way of seeing 
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these differences in experiments consists in measuring the 
degree of polarization in the schemes with and without 
the dichroic beam sphtter. In both cases the degree of 
polarization is related to the corresponding Schmidt pa- 
rameter: 



p(4); 



1 



Kpoi 



1, P? 



2 qb 



1 



K- 



2qb 



: 1. 

(83) 



Specifically, in two simple cases of subsections E and F, 
section V, Eqs. ([79|) - (|83|) give the following results. 

1) Ci = C4 = with + = 1, which corre- 

sponds to the original ququart's wave function 



B. 



V2 



(4) 
HV 



B. 



V2 



(4) 
VH- 



(84) 



\Bl- 



Bl\.P, 



2 qb 



l-\Bi- B 
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In this case C2 qb 

whereas KP°^ = 2, C^"' ^ \ \d+ 
and the dependence of all these parameters on \B^\ is 
shown in Fig. |4] for real values of B+ and The 



B+P - andPW = 0, 



1.5 



0.5 



2qb 




0.2 0.4 0.1 



\B. 



2) C4 = B+ = with |Ci|2 + |B_|2 = 1. The corre- 
sponding original ququart's wave function is given by Eq. 
(|69l) . The Schmidt parameter, concurrence and degree 
of polarization found in the two-qubit model are equal 

^ C2qb = and P2qb = 



now to K2qb = 



yet KP°^ 



2~\B., 

whereas in the case of mixed states we 
2 



p(4) 



2-2|B_|2 + |B_| 
1 — l-B-p. In this case, again, the concurrences C^"' and 
C2 qb coincide with each other and again because of a spe- 
cific choice of the ququart's constants C4 = B^ = 0. But 
the degree of polarization P^^^ does not coincide with 
P2 qb and the Schmidt parameter KP°^ does not coincide 
with K2qb- Dependence of all these quantities on the 
mixing parameter \B^\ is shown in Fig. [5l The dif- 



2 

1.5 
1 

0.5 


0.2 0.4 0.6 1 
Figure 5: The same as in Fig. |4]but for C4 = = 0. 

ference between the curves for the Schmidt parameter 
K{\B-\) and the degree of polarization P{\B^\) in cases 
of mixed state and of the two-qubit model is rather well 
pronounced and experimentally measurable. 













K2qb 















Figure 4: The Schmidt parameter KP°^ of the mixed state ll63t 
and the degree of polarization p(4) Qf ^^jjg state ift^) ||g4[| (sohd 
lines) and the same quantities for the two-qubit state I I74II (dashed 
Hnes) with coinciding curves for concurrencies (7''°' and C2 qb (solid 
line) for the case Ci = C4 = 0. 

differences between the curves corresponding to the 
ququart's mixed state and to the two qubit state are 
rather well pronounced in the case of the Schmidt 
parameter and the degree of polarization whereas the 
curves for the concurrences C^°' and C256 coincide. 
This last coincidence is accidental and, under the same 
conditions Ci = C'4 = 0, occurs only owing to the 
assumption that both constants B^ and B_ are real. If, 
for example, the constant B^ is complex rather than 
real, B^ = e^^jS+l (with real B-), then the concurrence 
of the mixed polarization state (j8ip remains the same 
as shown in Fig. (g]), C^"' = |l - 2|S_|2|, whereas the 
two-qubit concurrence (|79p becomes equal to C2qb = 
[(l-|i?„|2)2-2(l-|i?_|2)|i3_|2cos2^+|B_|4]^/^ 
and CP°' ^ C2 qb if only ^ ^ or tt. 



VII. CONCLUSION 

In Conclusion, biphoton polarization-frequency 
ququarts are pure states of two photons with two de- 
grees of freedom: polarization and frequency. Owing to 
indistinguishability of photons, their total wave functions 
must be symmetric with respect to the particle trans- 
positions. This gives rise to the intrinsic entanglement 
of such states related to their symmetry (symmetry 
entanglement). In addition, entanglement can arise 
also if the ququarts are formed by a superposition of 
several occupation-number configurations (configuration 
entanglement). Both symmetry and configuration 
entanglement are summed in the total entanglement of 
ququarts characterized by their Schmidt decomposition 
and Schmidt parameter K. If, however, either frequency- 
or polarization- variables are considered as "invisible" , 
the full density matrix of ququarts has to be averaged 
over these invisible variables. As a result one gets mixed 
two-qubit states depending only on either polarization 
or frequency variables. Such mixed states have both 
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classical and quantum correlations, and parameters 
characterizing levels of these correlations are found 
and analyzed. In the case of mixed states the Schmidt 
parameter K is shown to be not related to the degree 
of entanglement. But it remains useful as characterizing 
the degree of polarization of ququarts. The degree of 
polarization of ququarts is shown to be identical to the 
degree of polarization occurring in mixed polarization 
states arising after averaging of the full quqart's density 
matrix over the frequency variables. This degree of 
polarization can be found experimentally by means of 
a direct measurement of the Stokes parameters of the 
original ququart biphoton beam not split for two parts 
by DBS. On the other hand, in experiments with DBS, 
by measuring parameters of only one part of a split 
ququart 's beam, one can find the degree of polarization 
of such half a beam, different from the degree of polar- 



ization of the whole ququart 's beam. Such measurement 
permit to see two-qubit features of ququarts but they do 
not determine the degree of polarization of the ququart 
beam as a whole. In two simple cases the difference 
between measurements to be done with and without 
using the DBS is illustrated by pictures of Figs. |4]and[5] 
which show that the difference is rather well pronounced 
and is accessible for experimental control. 
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